UNIFORMITY OF HARMONIC MAP HEAT FLOW AT 
INFINITE TIME 



LONGZHI LIN 



Abstract. We show an energy convexity along any harmonic map heat flow 
with small initial energy and fixed boundary data on the unit 2-disk. In 
particular, this gives an affirmative answer to a question raised by W. Minicozzi 
asking whether such harmonic map heat flow converges uniformly in time 
strongly in the W /1 ' 2 -topology, as time goes to infinity, to the unique limiting 



harmonic map. 



1. Introduction 

Given a compact Ricmannian manifold AA and a closed (i.e., compact and with- 
out boundary) Ricmannian manifold Af which is an isometrically embedded sub- 
manifold of W 1 , we can define the Dirichlet energy of a map u £ W l ' 2 (A4, Af): 



(1.1) Energy(u) = E(u) = \ [ \Vu\ 2 dv M , 

where W 1 ' 2 (A4,Af) is the class of maps 

£ L} oc (M,M. n ) : / \Vu\ 2 dv M < +oo, u(x) £ Af a.e. x £ M 

JM 



>M 

The tension field t(u) £ T(u* (TAf)) is the vector field along u representing the 
negative L 2 -gradient of E(u) . A weakly harmonic map u from A4 to Af is a critical 
point of the energy functional E(u) in the distribution sense, i.e., the tension field 
t(u) vanishes, and it solves the Euler-Lagrange equation 

(1.2) - A M u = A(u)(Vu, Vu) , 

where u = (it 1 , ...,u n ) and A(u) denotes the second fundamental form of Af K™ 
at the point u. We refer to this system of elliptic equations as the harmonic map 
equation. 

A natural way to control the tension field for an energy minimizing sequence of 
maps and to get the existence of harmonic maps from AA to Af is to consider the 
initial(-boundary) value problem: 

{ut- A M u = A(u)(V«, Vit), onMx(0,T) 
u(x, 0) = ito(ar), for x £ M 

u(x,t) = x(x) = u \om, for all t > 0, x £ dM if dM ^ , 

where u = (u , u") and T > 0. We refer to this system of parabolic equations as 
the harmonic map heat flow, to the map uq as the initial data, and to the map x 
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as the boundary data. Given u <E W 1,2 (M t Af) and \ = u \qm G W^ 2 (dAi,Af), 
wc define u G W 1 ^ 2 (Ai x [0,T],7V) to be the weak solution of (1.3) if 

(1.4) f [ (u t ,0 + <V«,VO-<>4(«)(Vu,Vtt),Ottedt = 

JO JJK 

for any £ G C C °°(X x (0,T),R n ). 

In the fundamental paper [ES64] where the harmonic map heat flow was first 
introduced, Eells and Sampson proved that the harmonic map heat flow exists for 
all time in the case that the source domain Ai (of arbitrary dimensions) is without 
boundary and the target manifold Af has non-positive sectional curvature. They 
also proved that there exists some sequence of times U /*■ +oo such that 

Moo = hm u(-,U) 

i— >oo 

is a harmonic map from Ai to N '. The case in which the source domain Ai has 
boundary was dealt with by Hamilton in [Ha75] under the same curvature assump- 
tion on Af. The question of uniformity of the convergence in time of the flow 
considered by Eells and Sampson in [ES64] was left open at that stage, but it was 
settled later by Hartman in [Har67]. We shall state their results in the following 
theorem. 

Theorem 1.1 ([ES64, Har67]). Suppose that Ai and Af are two closed Riemann- 
ian manifolds and that Af has non-positive sectional curvature. Then given any 
uq G C l (Ai,Af), the harmonic map heat flow has a unique solution u £ C 1 (Ai x 
[0,oo),7V)nC°°(X x (0,oo),7V). Moreover, 

(1.5) = lim u(-,t) 

t— >oo 

exists uniformly in C k -topology for all k > and Uqo is a harmonic map homotopic 
to Uq. 

Other similar uniformity results were obtainable under various of assumptions on 
the target manifold Af, such as Af is real analytic (see Simon [Si83]) or Af admits a 
strictly convex function, see also Topping's interesting work [To97] in this direction 
for harmonic map heat flow in a special case in which both the source and target 
manifolds arc 2-sphcrcs S 2 . 

When the dimension of the source domain Ai is two it is particularly interest- 
ing because the energy functional E(u) and the harmonic map equation (1.2) arc 
conformally invariant in this critical dimension. Regarding the harmonic map heat 
flow (1.3) from surfaces to a general closed target manifold Af, the first fundamental 
work was for the case dAi = which was due to Struwe [Str85], where "bubbles" 
may occur and have been analyzed in detail. This result was then extended to the 
case dAi ^ with Dirichlet boundary condition by Chang in [Ch89]. If the initial 
energy E(uq) is sufficiently small, it is well-known by now that the weak solution 
of (1.3) is smooth (in the interior) by the results of Frcire [Fr95, Fr96] using the 
so called moving frame technique introduced by Helein (see e.g. [He02]). We state 
their e-regularity theorem that we shall require in this paper in the following. For 
the self-containcdncss of the paper, we will include an alternative proof of this e- 
regularity theorem using the main tool of our current work that we call Riviere's 
gauge decomposition, see Theorem 3.7. 
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Theorem 1.2 ([Fr95, Fr96], cf. [Str85, Ch89, Wall]). Let M be a simply-connected 
compact Riemannian surface andAf be a closed Riemannian manifold. There exists 
Eq > depending only on Ai and Af such that the following is true. For each initial 
data uo € W 1,2 (.M, Af) with E(uq) < Eq and the boundary data x = uq\qm in 
the case that dAi ^ 0, there exists a unique global weak solution u € W 1,2 (Ai x 
[0,oo),Af) for which E(u(-,t)) is non-increasing in t. Also, u is smooth in A4 x 
[1, oo ) and for any ti > t\ > 1 we have 

(1.6) 2 f 2 f M 2 =/ \Vu{-M)\ 2 ~l |V U (-,t 2 )| 2 . 

J t\ JM JBi JBi 

Moreover, there exists some sequence of times ti /* +oo such that 

(1.7) Moo = lim u(-,U) 

i—too 

exists in the C k -topology for any k > and Uoo is a harmonic map from Ai to Af. 

Remark 1.3. In particular, in order to avoid the "bubble" (singularity) along the 
harmonic map heat flow, a priori we may choose eo < Ki + K2 where 

K x = inf {E(v) \v e W h2 {A4,Af) and v\ dM = X } 

and 

K<i = inf {E(v) \v: § 2 — > Af is non-constant and harmonic} > . 

Remark 1.4. Frcire's regularity results for harmonic map heat flow in [Fr95, Fr96] 
is a parabolic version of the regularity theorem of Helein which states that weakly 
harmonic maps from surfaces are regular (see e.g. [He02]). 

A tempting question to ask is that, for a general closed target manifold Af (with- 
out imposing additional geometric assumptions on Af), whether one could establish 
uniformity results for the harmonic map heat flow similar to Theorem 1.1. In par- 
ticular, is the convergence (1.7) in Theorem 1.2 uniform for all time in the natural 
T4 /1,2 -topology, say? In view of the conformal-invariance of the energy functional 
E(u) in dimension two, the condition of small energy seems to be a natural candi- 
date to work with in order to get such uniformity of the convergence in time for the 
flow. We will show in the following that this is indeed the case. In what follows we 
will concentrate on the case that source domain Ai is a simply-connected compact 
Riemannian surface with boundary. More precisely we focus on domains which are 
conformally equivalent to the unit 2-disk B\ CK 2 . From now on we will only work 
on B\\ 

flit - Au = A(ii)(Vu, Vu), on B x x (0, T) 
u(x,0) = u (x), for x e Bi 

u(x, t) = x( x ) = uo\ob 1 , for all t > and x € dB\ , 

where A is the usual Laplacian A = Yli=i jfc? m • ^ ai 'g umen t s could be 
easily modified to apply to the general case. 

Notations 1.5. In what follows, V = (d x ,d y ) is the gradient operator in K 2 and 
V x = (—d y , d x ) denotes the orthogonal gradient (i.e., V -1 is the V-operator rotated 
by tt/2). 

Now we state the main theorem of this paper. 
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Theorem 1.6. Let Af be a closed Riemannian manifold. There exist so, Tq > 
depending only on Af such that if u £ W 1,2 (B± x [0, oo), Af) is a global weak solution 
to the harmonic map heat flow (1.8) with E(uq) < Eq, E(u(-,t)) is non-increasing 
in t and u(-,t)\gs 1 = X f or a ^ t > 0, then for all t2 > t\ > Tq we have the energy 
convexity 

(1.9) - f \Vu{-M)-Vu{-,t 2 )\ 2 < f \Vu{-M)\ 2 -f \Vu(-,t 2 )\ 2 . 

J B± J B\ J Bi 

Remark 1.7. We do not know if the energy convexity (1.9) holds for all i 2 > ii > 0, 
cf. [Wall]. In the following arguments we agree to let Eq be sufficiently small and 
To be sufficiently large as needed. 

Our approach to the proof of Theorem 1.6 is based on the technique that we call 
Riviere's gauge decomposition introduced by Riviere in [Ri07], see also Section 3. 
Immediate applications of Theorem 1.6 are the following corollaries. 

Corollary 1.8. Let Af be a closed Riemannian manifold. There exists Eq > 
depending only on Af such that if u £ W 1,2 (B± x [0, oo), Af) is a global weak solution 
to the harmonic map heat flow (1-8) with E(uq) < £q, E(u(-,t)) is non-increasing 
in t and u(-, t)\dBt = X f or a ^ t > 0, then 

(1.10) u{',t) — > "oo uniformly as t — > +oo strongly in W ' \B\, R") , 

where Uoo is the unique harmonic map with E(uoo) < Eq and boundary data X- 

Corollary 1.9. Let AA be a two dimensional domain that is conformally equiva- 
lent to Bi and has smooth boundary, and let Af be a closed Riemannian manifold. 
Suppose the initial energy E(uq) < Eq, then the harmonic map heat flow (1.3) with 
initial data uq € C 2a (A4, Af) and boundary data x € C 2a {dAA, Af), considered 
by Chang in [Ch89], converges uniformly in time strongly in W 1,2 {AA, Af) to the 
unique harmonic map € C^ a (A4,Af). 

Remark 1.10. We do not know if a harmonic map heat flow can be non-uniform 
without the small energy assumption. In view of the non-uniqueness results of 
Brezis and Coron [BC83] and Jost [Jo84] for harmonic maps (with large energy) 
sharing the same boundary data on dB\, it is quite possible that the small energy 
assumption is necessary for the energy convexity and uniform convergence of the 
flow in Theorem 1.6 and Corollary 1.8 to hold. 

Remark 1.11. In [CM08-1] Colding and Minicozzi showed an energy convexity for 
weakly harmonic maps with small energy on B\: there exists eo > such that if 
u,v G W ' (Bi,Af) with u\qb 1 = v\dBn E(u) < £o, and u is weakly harmonic, then 
we have the energy convexity 



(1.11) 




See our recent work [LL12] for an alternative proof of this energy convexity using 
the same techniques used in the present paper. A direct consequence of (1.11) is 
that Uoo in Corollary 1.8 is unique in the class 

{v € W 1 ' 2 (Bi,R n ) : E(v) < e and v\ dBl = X } , 

see [CM08-1, Corollary 3.3]. 
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The paper is organized as follows. In Section 2 we present some heuristic argu- 
ments and elaborate on the idea of the proof of the main Theorem 1.6. In Section 
3 we review the main tool of our proof, namely, Riviere's gauge decomposition 
technique adapted to the case of harmonic map heat flow 1 . In Section 4 we show 
improved estimates for Riviere's matrices B and P, which are the two key ingredi- 
ents of our proof. We finish the proof of our main theorem in Section 5. 

Acknowledgement. The author would like to thank Professor William Minicozzi 
for his continued guidance and support during the author's Ph.D. study at The 
Johns Hopkins University. The author would also like to thank his collaborator 
Tobias Lamm for all the valuable discussions. 



2. Heuristic arguments and the idea of the proof 



In this section we will present some heuristic arguments and sketch the basic 
idea of the proof of Theorem 1.6. We will abbreviate u(-,t) to u(t). In order to 
prove the energy convexity (1.9) along the harmonic map heat flow, i.e., there exists 
some To > such that for all t 2 > h > To we have 

(2.1) -j \Vu{-M)-Vu{-M)\ 2 < f |Vn(-,*x)| 2 - f \Vu(-,t 2 )\\ 

J B\ J B± J B\ 

it suffices to show 

(2.2) * > - (7 |V^i)| 2 - / |V^ 2 )| 2 ) -if \Vu(h) Vu(t 2 )\ 2 , 

\J B\ J B\ / J B\ 

where (using that u{-,t)\aB-i = x f° r au t >Q and the flow equation (1.8)) 
*:=/ IV^)! 2 -/ \Vu{t 2 )\ 2 -f \Vu(h)-Vu(t 2 )\ 2 

J B\ J B\ j B\ 



(V«(ti)-Vu(t2),V«(t2)) 
(2.3) = - 2 / (u(tx) - u(f a ), u t {t 2 ) - A(«)(V«, Vu)(t 2 )) • 



Now note that for any p,q G J\f, there exists some constant C > depending only 
on M such that | (p — q) 1 - 1 < C\p — q\ 2 , where the superscript _L denotes the normal 
component of a vector (see e.g. [CM08-2, Lemma A.l]). Therefore, using the fact 
that A(u)(Vu, Vu) _L T U M and Cauchy-Schwarz inequality, (2.3) yields 



tf>-2( \u(h) - u(t 2 )\' j ^ \u t (t 2 )\ 2 j -CJ^ \{u{tr)-u{t 2 )) x \\Vu{h)\< 



>-2Vh=h( / \u t \ 2 \ / \u t (t 2 )\ 2 \ -C Hh)- u(t 2 )\ 2 \Vu(t 2 )\ 

\Jti JB X J \J Bx J JBx 

where we also used the smoothness and compactness of the target manifold M . 
Here and throughout the rest of the paper, C > will denote a universal constant 
depending only on Af unless otherwise stated. 



^Alternatively, one may use the results of Freire [Fr95, Fr96] where Helein's estimates involving 
moving frame were used, however, Riviere's approach seems to be more elementary and self- 
contained. 
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Since we have (1.6) and £o can always be chosen sufficiently small, we know that 
(2.2) will be achieved if we can show the following two key propositions. 

Proposition 2.1. Let u(x,t) be as in Theorem 1.6, then there exists Tq > such 
that for all t 2 > t\ > To we have 

I I , , 2 



(2-4) / |«*C*a)r < T T / / ^ 

Remark 2.2. The key point of Proposition 2.1 is that (2.4) is valid for all t 2 > t\ > 
T . We will see that, in fact, J B |u t (t)| 2 is non-increasing along the flow after T , 
which yields (2.4), cf. Lemma 2.5 and (5.16) below. A similar but weaker estimate 
was shown in [Str85] when the source domain of the heat flow is boundarylcss 
(see equation (3.5) of [Str85]), which turned out to be the key estimate needed in 
Struwe's proof. 

Proposition 2.3. Let u{x,t) be as in Theorem 1.6, then there exists Tq > such 
that for all t 2 > t\ > To we have 



(2.5) / \u{t{)-u{t 2 )\ 2 \Vu{t 2 )\ 2 < Cs 



[ |V«(*i) - V«(ta) 

JBi 



If one was able to get 
(2.6) \\Vu(t 2 )\\ L oo (Bl) <C^, 

then (2.5) would have been automatically true by Poincare's inequality. However, 
without imposing any regularity information on the boundary data x-i it will be 
hopeless to get such strong global pointwisc gradient estimate. In fact, even if we 
look at the stationary case, i.e., W 1,2 -weakly harmonic maps on Bi, it is easy to 
convince oneself that it is unreasonable to expect regularity with global estimates 
on the whole B\ better than W 2,2 in general. 

Nevertheless, not all hope is lost to show estimates (2.4) and (2.5). Indeed, the 
following lemma is true which validates Proposition 2.3 under some extra assump- 
tion. 



Lemma 2.4. Let u(x,t) be as in Theorem 1.6 and suppose that for all t 2 > ti > 

rl, 

Aip = \Vu{t 2 )\ 2 in B 1 , 
tp = on dBi , 



To > 1 we can solve the following Dirichlet problem for ip £ W ' PI L°°(Bi): 
(2.7) 



with estimate 

(2.8) IMU-(Bx) + IIV^Mbo < Csq. 
Then Proposition 2.3 holds. 

Proof. The proof is essentially taken from [CM08-1]. Substituting (2.7) into the 
left-hand side of (2.5) yields (using also that u(ti) = u(t 2 ) = \ on dB\) 

\u{h) - u(t 2 )\ 2 \Vu(t 2 )\ 2 = [ M<i)-u(t2)| 2 A</>< / |V|u(ii)-w(t 2 )| 2 ||VV| 

JB-l JB-l 

(2.9) <2(7 |Vu(fi) - Vu(t 2 )A ( [ Kti)- W (t 2 )| 2 |WfV , 
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where we have applied Stokes' theorem to div(\u(ti) — u{t 2 )\ 2 V ip) and used Cauchy- 
Schwarz inequality. Now applying Stokes' theorem to div{\u{ti)-u(t 2 )\ 2 ipV4>) and 
using that Ai/j > and (2.9), we have 

|u(ti)-u(f 2 )| 2 |V^| 2 < f MHh) - u(t 2 )\ 2 A^ + |V|«(ti) - u(t 2 )\ 2 \\V^\) 

(2.10) <4||V||ioc Qf \Vu(h)-\7u(t 2 )\ 2 ^ ' ( [ \u(h)-u(t 2 )\ 2 \Vi;\ ^ ' 
and so that 



Bi 



(2.11) |u(ti)-«(* 2 )r|V^| 2 ^ 7 < M\i/>\\l- \Vu(h)-Vu(t 2 )\ 2X ' 

Finally, substituting (2.11) back into (2.9) and combining with (2.8) (and choosing 
£o sufficiently small) yield 

\u(h) - u(t 2 )\ 2 \Wu(t 2 )\ 2 < CU\\ L oo f \Vu(h) - Vu(t 2 )\ 2 

J Bi 

<Ce a [ |Vu(ti)- Vu{t 2 )\ 2 , 

J B\ 

which is just (2.5). □ 

Similarly we can show the following lemma which states, under some extra con- 
dition, J B \u t (t)\ 2 is non-increasing along the harmonic map heat flow after some 
To > and Proposition 2.1 could be validated in this case. 

Lemma 2.5. Let u(x,t) be as in Theorem 1.6. For any t 2 > t\ > To > 1, 
suppose that for any to € [ti,t 2 ] we can solve the following Dirichlet problem for 
-0 € W 1,2 n L°°(Bi): 

(Aip = \\/u(t )\ 2 m Si, 
1 tp = on dB!, 



(2.12) 



wif/i estimate 

(2.13) IMU~( Bl ) + IIVV'II^bo < Ce . 
TTien we /iave 

(2-14) / \u t (t 2 )\ 2 < f \u t ( tl )\ 2 . 

J Bi JBi 

In particular, Proposition 2.1 holds if (2.12) and (2.13) are valid for any to € [ti,i2] 
and any t 2 > t\ > To > 1 . 

Proof. Differentiate the flow equation (1.8) with respect to t, multiply with u t , 
and integrate over B\ x [£1,^2], we have (e.g. treating u t as a difference quotient: 
u t (-,t) = lim /l ^ + i u {'it + h) — u(-,t))/h which is zero on dB\ for all t > 1) 

Iff d t \u t \ 2 + f 2 / |Vu t | 2 < C f f \u t \ 2 \Vu\ 2 + \u t \\Vu\\Vu t \ 

Z JU JBi Jti JB, Ju JBi 



• t\ J B± J t\ J B\ J t\ J B 

r t 2 r rt 



(2-15) <\l 



ti JB 
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Since (2.12) and (2.13) are valid for any to G [ti, £2], we can use the same argument 
as in the proof of Lemma 2.4 to get estimate for 

/ M 2 |Vu| 2 

at the time to slice. Indeed, similar to (2.5) (i.e., replacing u(ti) — ufa) by Ut(to)), 
for any to G [ti, £2] we have 

(2.16) f \u t \ 2 \Vu\ 2 (t ) < CU\\ L oo [ \Vu t (to)\ 2 <Ce f \^u t {t )\ 2 . 

J B\ J B\ J B\ 

Inserting (2.16) back into (2.15) (for any to 6 [ti,ta]) we see that the right-hand 
side of (2.15) can be absorbed into the the left-hand side if we choose £0 sufficiently 
small. This implies that we have (2.14) for any such ti > t\ > To. 

If (2.12) and (2.13) arc valid for any to G [h,t 2 ] and any t 2 > h > T > 1, 
then, in view of (2.14), estimating by the mean value of \ut\ 2 over B\ x [ti, t?\ gives 
Proposition 2.1. □ 



Therefore, everything boils down to validating the assumptions in Lemmas 2.4 
and 2.5, i.e., the existence of such functions ip's satisfying (2.7), (2.8) and (2.12), 
(2.13), respectively, for any t > T for some T > 1. We point out that, a prior 
we only know that the energy density |Vu(£)| 2 lies in L 1 (Bi) with global estimate 
|| |V«(f)| 2 < £ o f° r an y fixed t. But L 1 is the borderline case in which the 
standard L p -theory for the Dirichlet problem (2.7) with estimate (2.8) fails! 

However, the following regularity theorem for boundary value problems in the 
local Hardy space /i 1 (i?i) sheds a new light to validate the assumptions in Lemmas 
2.4 and 2.5. Here the local Hardy space /i 1 (i?i) is a strict subspace of L 1 (Bi) and 
we will recall its definition in Definition 2.8 below. 

Theorem 2.6. (cf. [Sc94, Theorem 1.100] and [CKS93, Theorem 5.1}) Let f G 
/i 1 (i?i) and assume that f > a.e. in B\. Then there exists a function tp G 
i°° f~l W Q ' (B±) solving the Dirichlet problem 

( Aip = f in B x , 
(2 - 17) ( 4, = on &Bi . 

Moreover, there exists a constant C > such that 

(2-18) IMU~ (Bl) + IIVV'II^bo < C ||/|Ui( Bl ) . 



Proof. An elementary proof can be found in [LL12, Theorem A. 4], which follows 
along the lines of a corresponding result in [Sc94]. □ 

Remark 2.7. This theorem can be thought of as a generalization of the Wente's 
lemma 3.6 in Section 3, cf. [Mu90]. For a more general version of this theorem we 
refer to Chang, Krantz and Stein's work in [CKS93]. 

Definition 2.8. ([Mi90]) Choose a Schwartz function </> G Qf(Bi) such that 
J B 4>dx = 1 and let (j>t(x) = t _2 0(j). For a measurable function / defined in 
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B\ we say that / lies in the local Hardy space h}{B\) if the radial maximal func- 
tion of / 

(a;) = sup \<k*f\ ( x ) 

0<t<l-\x\ 

belongs to L 1 (i?i) and we define 

(2-20) H/lkKBi) = Wf*(x)\\ LHBl y 

It follows immediately that is a strict subspacc of L l {B\) and H/IIli^) < 

ll/IU 1 (Bi)- It is also clear that if / e L p (Bi) for some p > 1 then H/H;,.!^) < 

C||/IUp(bi)- 

We shall remark that the local Hardy spaces h 1 (or the global version H 1 ) act 
as replacements to L 1 in Caldcron-Zygmund estimates. Therefore by Theorem 
2.6, if we can somehow manage to obtain a "slightly" improved global estimate 
for |Vw| 2 from L l {Bi) to for all to > T Q , it will be sufficient to validate 

the assumptions in Lemmas 2.4 and 2.5. As mentioned above, the subtlety is 
that, without imposing any regularity information on the boundary data \i global 
estimates are very difficult to obtain. 

The rest of the paper is devoted to validating the assumptions in Lemmas 2.4 
and 2.5. Namely, in view of Theorem 2.6, it suffices to show there exists To > 
such that 

(2.21) |||Vii(t )| 2 ||/ l i(i3 1 ) < Ceo for any t > T . 

The point here is that no pointwise estimate on Vw such as (2.6) is needed, and 
instead, a (weaker) improved global integral estimate (2.21) will be sufficient and 
turns out to be the key to the proof of the main Theorem 1.6. 

3. Analysis of harmonic map heat flow using Riviere's gauge 

Regarding the regularity of weakly harmonic maps from surfaces, Helein (see 
e.g. [He02]) proved the interior regularity with the help of the so called Coulomb 
or moving frame, and Qing [Qi95] showed the continuity up to the boundary in 
the case of continuous boundary data based on Helein's technique. More recently, 
Riviere [Ri07] succeeded in writing the 2-dimensional conformally invariant non- 
linear system of elliptic PDE's (which includes the weakly harmonic map equation 
(1.2)) in the following form: 

(3.1) -Au i =n)-Vu j i=l,2,...,n or -Au = fi-Vu, 

with fi = (fij-)i<*,i<n G L 2 (B 1 ,so(n) g) A 1 ]^ 2 ) and fK = -0? (antisymmetry). 
Here and throughout the paper, the Einstein summation convention is used. We 
refer to the system of equations (3.1) as Riviere's equation. In particular, this 
special form of the nonlinearity enabled Riviere to obtain a conservation law for 
this system of PDE's (see (3.8) below), which is accomplished via a technique that 
we call Riviere's gauge decomposition. More precisely, following the strategy of 
Uhlenbeck in [Uh82], Riviere [Ri07] used the algebraic feature of f2, namely f2 being 
antisymmetric, to construct £ € W l,2 {Bi 1 so(n)) and a gauge transformation matrix 
P E W 1,2 n L°°(Bi, SO(n)) (which pointwise almost everywhere is an orthogonal 
matrix in K™ xn ) satisfying some good properties. 



(2.19) f*{x) 



sup 

0<i<l-|x| 



B t {x) 



1 



x-y 



f(y)dy 
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Theorem 3.1. (|Ri07, Lemma A.3]) There exist e > and C > such that for 
every ft in L 2 (Pi, so{n) ® A 1 M 2 ) satisfying 

! m 2 < e, 

there exist £ g W^ 1 ' 2 (Pi, so(n)) and P g W 1 ' 2 (B 1 , SO(n)) such that 

(3.2) = P T VP + P T flP inB 1 with £ = on dB t , 
and 

(3.3) l|Vf||L»( Bl ) + IIVP|U= (Bl ) < c\\n\\ LHBl) . 

Here the superscript T denotes the transpose of a matrix. 

Remark 3.2. Multiplying both sides of equation (3.2) by P from the left gives that 
(with indices and 1 < m, z < n) 

(3.4) VP] = P,; V-%" 1 - Qi P/, 1 < i, j < n . 

Remark 3.3. Besides Uhlcnbeck's method there is another way to construct the 
gauge tranformation matrix P, namely one can minimize the energy functional 

(3.5) E(R) = [ |P T VP + R T flR\ 2 
among all R g W ia {B l , SO{n)), see e.g. [Cho95] and [Sc09]. 

Another key result from Riviere's work is the following theorem, which was 
proved based on Theorem 3.1. 

Theorem 3.4. ([Ri07, Theorem IA\) There exist e > and C > such that for 
every £1 in L 2 (B\, so(n) ® A 1 ^) satisfying 

{ |ft| 2 < e, 

JBi 

there exist A g W 1 - 2 n C°(Bi,Gl n (M)), A = {A+Id) P T g L°° n W 1 ' 2 ^!, GZ„(R)) 
and P g W / 1 ' 2 (Pi,M„(E)) sucft ffta* 

(3.6) VA- Afl = V^P 
and 

(3.7) \\A\\ w ^ {Bl) + \\A\\ L ~ (Bl) + \\B\\ wl , 2{Bl} < C||0|| L2(fll) . 

Remark 3.5. Combining (3.6) with (3.1) one obtains the conservation law (in the 
distribution sense) for Riviere's equation (3.1): 

(3.8) div (AVu + BV^u) = . 

Equation (3.1), first considered in such generality by Riviere in [Ri07], generalizes 
a number of interesting equations appearing naturally in geometry, including the 
harmonic map equation (1.2)), the P-surfacc equation and, more generally, the 
Eulcr-Lagrangc equation of any conformally invariant elliptic Lagrangian which is 
quadratic in the gradient. We remark that the harmonic map equation (1.2) can 
be written in the form of (3.1) if we set 

(3.9) := (n*-)i<ij<n where Oj := [A\u) hl - A 3 '(tt)«,i]W . 
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A central issue is the regularity of the weak solution u to this system of equations 
(3.1). Based on the conservation law (3.8), Riviere proved the (interior) continuity 
of any W 1,2 weak solution u to (3.1). This also resolved two conjectures by Heinz 
and Hildcbrandt respectively, see [Ri07] . Wc point out that the harmonic map heat 
flow (1.8) on B\ can be written in the form 

(3.10) u t - Am = fi • Vit on Bi x (0, T) , 

where Vt is as in (3.9). 

The deep reason for Riviere's argument to work is that once the conservation 
law (3.8) is established, then equation (3.1) can be rewritten in the form 

div(AVit) = V^B • Vit, 

The right hand side of this new equation lies in the Hardy space H 1 by a result 
of Coifman, Lions, Meyer and Semmes [CLMS93]. Moreover, using a Hodge de- 
composition argument, one can show that u lies locally in W 2,1 which embeds into 
C° in two dimensions, cf. the proof of Theorem 3.7 below. The key to this fact 
is a special "compensation phenomena" for Jacobian determinants, which was first 
observed by Wentc in [We69]. We will refer it to the following Wente's lemma, for 
which an elementary proof can be found in [BC83] and [Hc02, Theorem 3.1.2]. It 
will also be the key ingredient of our proof in this paper. 

Lemma 3.6 (Wente's lemma [Wc69]). If a, b € W 1 ' 2 (Bi,R) and w be the solution 
of 

da db da db 



(3.11) 



oy ox ox oy 
dw 

or — — = on dB\ 

ov 



Then w £ C° n W ' (Bi,M.) and the following estimate holds 

(3.12) IM| £ o. (Bl) + HVH|i2 (Sl) < C , ||Va|| L 2 (Bl) ||V6|| L 2 (Bl) , 
where we choose J g w = for the Neumann boundary data. 

Now let u(x,t) € W 1 ' 2 (B 1 x [0, oo),7V) be a global weak solution to the har- 
monic map heat flow (1.8) with E(uo) < £oi E(u(-,t)) is non-increasing in t and 
""(•7 0l<9Bi = X f° r all <: > as in Theorem 1.6. First note that for a.e. to G (0, oo) 
we have Ut(to) e L 2 (Bi). Then for any fixed to such that u t (to) E L 2 (Bi), as in 
(3.9), we have 

O(to) = (O}(t ))i<i,j< n where Slfa) = [A l (u(t ))jj - (u(to)) t j]Vu l (t ). 
We will abbreviate Q(t ) to A(u(t ))V u(t ) . Moreover, 

(3.13) / \n(to)\ 2 < CE(u(t )) < Ceo. 

J Bi 

Therefore Riviere's Theorems 3.1 and 3.4 on the existence of gauge apply to this 
time to slice, and we find the existence of matrices P(to) G W 1,2 (Bi, SO(n)), 

A(to) = (A(t Q ) + Id)P T (to) G L°° n W 1,2 (Bi,Gl n (R)) 

and 

B(to) G W 1 ' 2 ( J Bi,M„(R)) 
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such that 

(3.14) VA(t ) - A(to)Q(to) = V x B(t ) 

with the corresponding estimates (3.3) and (3.7). 

Combining (3.14) with the harmonic map heat flow equation (3.10) yields (omit- 
ting the index to) 

div (AVu + BV^u) = VA • V« + AAu + VB ■ V ± u 

= VA ■ Vu + A(-Q • Vu + u t ) + VB ■ V ± u 

(3.15) = VA ■ Vu + (y ± B - VA) ■ Vu + Au t + VB • V x u 

= Au t . 

We refer to (3.15) as an "almost" conservation law. By the results of [CLMS93] and 
the standard L p theory, (3.15) readily implies that u(to) £ C°(Bi,M. n ). In fact, we 
have the following e-rcgularity theorem. 

Theorem 3.7. There exist eo > depending only onM such that if u £ W 1,2 (Bi x 
[0, oo),Af) is a global weak solution to the harmonic map heat flow (1.8) with 
E(uq) < £o, E(u(-,t)) is non-increasing in t and M(-,t)|as 1 = x f or a ^ i ^ 0, 
then u £ C°°{B 1 x [l,oo),7V). 

Proof. For any fixed to such that u t (to) £ L 2 (Bi), by Hodge decomposition (see 
e.g. [IM01, Corollary 10.5.1]), there exist D(t ),E(t ) £ W 1 - 2 {B 1 ,W l ) such that 
(omitting the index to) 

(3.16) AVu = VD + V ± E . 
Note that (3.15) implies 

( div (AVu) = -VB ■ V 1 ?! + Au t 
\ curl (AVu) = V ± A ■ Vu . 
Combining (3.16) and (3.17) we have 

f AD = -VB ■ V x u + Au t 
1 AE = V^A • Vu . 



(3.17) 



(3.18) 



Then by the results of [CLMS93] and via an extension argument, using the fact 
that Au t (t Q ) £ £ 2 (Bi), wc get AVu{t ) £ W}£{Bi). Therefore u(t ) £ 
which embeds into C°(Bi). 

Indeed, AVu(to) £ W,' (Bx) implies immediately that 

(3.19) Q(t ) = A(u(t ))V«(t ) G W£(Bi) . 

Then by [Ri09, Theorem III.4], the flow equation (3.10), (3.19) and the fact that 
u t (to) £ L 2 (Bi) yield that Vu(t ) £ Lf oc (Bi) for some p > 2. Note that this is 
valid for a.e. to £ (0, oo). Then via a standard bootstrapping argument we have 
Vu e Bj oc (B\ x [1,T]) for all q > 1 and any T > 1, see e.g. [Li96], and all higher 
order interior regularity follows. □ 

Again, we see that the "compensation phenomenon" enjoyed by the special Jacobian 
structure (see Wente's lemma 3.6) has played an important role here, and these 
Wente type estimates have many interesting applications as observed already in 
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e.g. [We69, BC83, BC84, Ta85, CLMS93, Hc02, Ri07, Ri08-1, Ri08-2] and our 
recent work [LL12]. 

4. Improved estimates on the matrices B and P 

Our main observation in this section are two hidden Jacobian structures for AB 
and AP, which are valid only for fi being in some special form, i.e., in our case we 
have n = A(u)Vu . This observation allows us to gain improved global estimate on 
the matrix B and improved local estimate on the matrix P. 

4.1. Improved global estimate on the matrix B. We will first show an im- 
proved global estimate on the matrix B. 

Proposition 4.1. Let u(x,t) be as in Theorem 1.6, then for any to 6 [l,oo) we 
have 

(4.1) ||fl(to)|| L o. (fll) < C f \Vu(t Q )\ 2 < Ce . 

Proof. We recall that fl is given by A[u)Vu as in (3.9) and therefore || r2(£o) II — 
Ceo for all to > 1. Now let £o be so small that Theorems 3.1 and 3.4 apply. Taking 
the curl on both sides of equation (3.14) yields 

(4.2) AB(t ) = -cnrl(A(t )A(u(t ))Vu{to)) = -Vu(t ) ■ V- L (A(i )A(u(i ))) . 

Combining the Jacobian structure of the right-hand side of (4.2) with the zero 
boundary condition of B and estimates (3.3) and (3.7), the Wente's lemma 3.6 
gives (4.1). Here we have also used E{u(t)) < e$ for all t > , the smoothness and 
compactness of the target manifold Af. □ 

4.2. Improved local estimate on the matrix P. We next show that AP also 
has a special Jacobian structure. 

Lemma 4.2. Let u(x,t) be as in Theorem 1.6, then for any to 6 [l,oo) such 
that u t {-,t ) G L 2 {B l ), there exist f(i ) G W 1 ' 2 {B 1 , so(n)), r){t ) € W 1 ' 2 (B 1 ,W l ), 
((t ) G W 2 ' 2 (Bi,]R' 1 ) and Q k (t ),Rk(t ) G W 1 ' 2 n L°°(B U Gl n {R)), k = l,...,n 
with 

||V£(to)|U2( Bl ) + ||Vf7(to)|Ua( Bl ) + ||VC(to) 11^(230 

+ (ll V( 9fc(*o)IU=(B 1 ) + W^RMWlhb,)) < CJTo 

k 

and 

(4-3) ||C(*o)lk^( Sl ) < C\\u t (t )\\ L 2 (Bl) 

such that 

AP(io) =VP(to) • V-'-f (*o) + VQ fc (to) • V V(*o) 
(4.4) + VR k (t Q ) ■ V ± u k (to) + div {Q k {to)VC k {to)) ■ 
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Proof. We omit the index to in the proof. Hodge decomposition and the estimates 
for the L°°-norms of A and B imply the existence of 77 £ W 1,2 {B\,W l ) and £ <E 
W 1,2 (Si,R") such that 

(4.5) V- 1 ?? + VC = AVu + BV 1 ^ 
with 

(4.6) IIVj/IU^bo + ||VC|U= (fll) < GIlVulU^) < CJT Q . 
Moreover, by the "almost" conservation law (3.15) we have 

AC = Au t e £ 2 (Bi) and Clsfli =0, 

which gives (4.3) by the standard L p -theory. Multiplying both sides of equation 
(4.5) by A^ 1 from the left gives (with indices) 

(4.7) W = (A~ 1 ) l k V ± r) k - (A- 1 B)' k V ± u k + (A" 1 )^* , 1 = 1,2, .... n . 
Taking divergence on both sides of equation (3.4) yields 

(4.8) AP] = VP4 ■ V^f - div (fi* P/) , 1 < i, i < n . 

Since = L4 4 (m) 2 . ; - ^(u) ij; ]Vu', combining (4.7) and (4.8) gives 
APj = VP^ • V^C™ - div [{A\u) z , - A*( M ) M ) ((A- 1 )^ V 

- (A _1 P)[.V ± u fe + ( J 4" 1 ) i fe VC fc ) P/] 

(4.9) = VP^ • V^f - V [(A>)^ - A*( U ) M )Pf (A" 1 )'] ■ V V 

+ V - A z ( M ),, i )P/(A- 1 P) / fe ] • V^ fc 

-div [((^( U )^-^( U ) M )Pf(^- 1 )i)VC fe ] ■ 

Defining (Q fc )j = -(^(u),,, -A'Cu)*,,).^- 1 ^ and (P fe )j = (A^u)^-A z (u). u y 
PHA~ 1 B) l k , where 1 < k,i,j < n, completes the proof. □ 



Next we prove a local estimate on the oscillation of the matrix P based on Lemma 
4.2. A key observation here is that whether a function is in the local Hardy space 
h}(B\) essentially depends on its local behavior (see Definition 2.8). This local 
oscillation estimate on P provides important information that we need to control 
the local behavior of |Vu| 2 . This point will become apparent in Section 5. As we 
shall see, the Jacobian structure of AP enters in a crucial way. 

Lemma 4.3. Let u{x,t) be as in Theorem 1.6, then for any to £ [l,oo) such that 
u t (-, to) G L 2 (Bi), any x G B\, any r > such that B2 r (x) C B\ and any y £ B r (x) 
we have 

(4.10) \P(y,t ) - P(x,t )\ < C(^+\\u t (to)\\ L 2 {Bl) ) . 

Proof. We will omit the index t in the proof. Let P € W 1 ' 2 (B 1 , M n (E)) be the 
weak solution of 

J AP = VP ■ V-^ + VQ fc ■ V-V + VPfe ■ V ± u k + div (Q fc VC fe ) in B l , 
\ P = on dB 1 , 

where Qk, Rk, r] k and ( k, s are from Lemma 4.2. 
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Then by Wente's lemma 3.6 and the standard L p -theory (and W 2 ' 2 (B 1 ) 
C°(Pi)), we have P G C°{B 1 ,M n {M.)) and 

(4-11) ||-P||z-(BO + l|VP||i2 (Bl) < C (e„ + ||ttt(to)IUa(Bo) ■ 

Since A(P-P) = in B t , wc know that V = P-P E C°°(Bi,M n (R)) is harmonic. 
Now for any igBi, any r > such that B2 r (x) C B\ and any y € B r {x) we have 

\V(y)-V(x)\ < Cr||W|| L » (Br(x)) < Cr^\\ VV\\ L , {B2r(x)) 

(4.12) < CUVFII^^^)) < C (||VP|| L2(B2r(x)) + ||VP|| ia(Bar(!e)) ) 

where we have used the mean value property of V and (4.11), (3.3). Combining 
(4.11) and (4.12) yields that for any x € Pi, any r > such that B2 r (x) C B\ and 
any y € B r (x) we have 

(4.13) \P(y,t ) - P(x,t )\ < C(^+\\ut(t )\\ L 2 (Bl) ) , 

which gives the desired estimate (4.10). □ 

5. VALIDATING (2.21) AND PROOF OF THEOREM 1.6 

With the results so far at our disposal, we are now in a position to validate (2.21). 
As mention above, the local estimate on the oscillation of the transformation matrix 
P in Lemma 4.3 will be the key ingredient. 

Lemma 5.1. Let u{x,t) be as in Theorem 1.6, then for any to € [l,oo) such that 
IK(*o)||z,*(Bi) < we have 

(5.1) \Vu{t )\ 2 e h 1 (B 1 ) 
with estimate 

(5.2) \\\Vu(t )\ 2 \\ hHBl) < Ce - 

Remark 5.2. Lemma 5.1 continues to hold if we replace f2*- = [A l (u)jj — A> (u)i^i\Vu l 
in (3.9) by a more general fl in the form fij = X)]Li f)i^ ul + ff}iV~ L tt', see [LL12]. 
Moreover, the condition ||itt(io)IU 2 (Bi) < £ o can be replaced by that ||ut(io)||.z>(Bi) 
is sufficiently small for some p > 1 . 

Proof, (of Lemma 5.1) By the assumption ||Mt(io)||i 2 (s 1 ) < \fto~ and Lemma 4.3, 
for any x € B\ , any r > such that i?2r (x) C B\ and any y & B r (x) we have 

(5.3) |P(^ )-P(Mo)| < C^o". 

We will omit the index to from now on. By Proposition 4.1, Theorem 3.1 and 
Theorem 3.4, for any x £ Si, any r > such that B2,(x) C Pi and any y € B r (x) 
we have (choosing eo sufficiently small) 

< ^|Vu| 2 (y) < (AVu + PV^u) ■ (P T \7u)(y) 

(5.4) = (AWu + BV^u) ■ [(P T (x) + (P T - P T (x))) Vu] (y) , 
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and therefore by (4.5) and (5.3) 

(V X 7? + VC) ■ (P T (x)Vu) (y) = (AVu + BV^u) ■ (P T (x)Vu) (y) 

(5.5) > \\Vu\\y) - {AVu + BV^u) ■ [{P T - P T (x)) Vu] (y) > \\Vu\ 2 (y) . 
Now we choose a function 

(5.6) <f> e C °°(Bi) with > 0, spt(0) C Bi , 4> = 2 on Ba, and / <j>dx=l. 

Moreover, we additionally assume that || V^Hl^^) < 100. Using (4.3) and (5.5), 
one verifies directly that (by Definition 2.8) 



|||Vu| 2 || h i (Bl) = / sup <j) t *\\7u\ 2 dx 

JB! 0<t<l— |je| 



<4 / sup 4> t * ({V^ri + VC) • (F T (i)Vu)) dx 

J Bi 0<t<l-|x| 

=4 / sup <j> t * [{P T {x)) l0 (V"V ■ S7u j + VC • Vm 7 )] dx 

JB 1 0<t<l-\x\ 
n 

(yv^^ ■ v^'|Ux (Sl) + ||vci| wl , 2(Sl )ll v ^IU 2 c^)) 

i,j=l 

<C||V- L 77|| L 2 (Bl) ||Vu|| L 2 (Bl) + Cy/e^\\u t \\ L 2 {Bl) < Ce , 
where we have used the facts 

(1) VV • W e h l {B x ) and ||VV ■ W\\ h i (Bl) < C\\Vr,\\ LHBl) \\Vu\\ L 2 iBl) 
for all i, j = 1, 2, n; 

(2) ||VC • W|| LP(Bl) < C||VC*||Hri, a(Bl) ||V«»'||L»(flO for any 1 < p < 2 and 
||/||fti(Bi) ^ C|l/IU"(Bi) for anyp > 1. 

To see (1), we first extend rf — rg^-r J B ^ rf and u J — rg^-r J Bj m j from B\ to R 2 which 
yields the existence of rf ,v? G W^' 2 (R 2 ) such that 

(5.7) f \Vf) l \ 2 <c( \Vrf\ 2 and f \VV\ 2 <c( \Vii?\ 2 

JM 2 JBi JR 2 J B\ 

and 

(5.8) Vrf = V?f and Vu J = Vit J a.e. in B a . 
Then by the results of [CLMS93] we know that 

|| • V^||„ 1(R2) : = / sup sup / 1^ {~~T^\ ( V± ^ ' VlV ) (») d W dx 

jR 2 <p£Tt>0 JBtix) 1 \ 1 J 

(5.9) < C||V?7l L 2 (E2) ||V^|| L 2 (R2) < C||V»?||L»(fl 1 )||V«|U a( fl I ) , 

where T = {</> € C°°(R 2 ) : spt(<£) C Si and ||V0||z,°° < 100}. By (5.8), (5.9) and 
Definition 2.8, it is clear that 

||VV-V^'|| fc i (fll) = ||V- L ^-Vu J '|| h i (fll) 
(5-10) < ||V^77 £ ■ V^||„i (a2) < GHVf/IU^BollVulUa^). 

This completes the proof of the lemma. □ 
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Now since u(x,t) G W 12 n C°°{Bi x [l,oo),7V) and the energy E(u(-,t)) is 
non-increasing along the flow as shown in (1.6), there exists T) > 1 such that 

(5.11) K(T )|| i3(fll) <V^. 

Then by Lemma 5.1 we know that |Vit(To)| 2 G /i 1 (Bi) with estimate 

(5.12) \\\Vu(T Q )\ 2 \\ hl{Bl) <Ce . 

Therefore, in view of Lemma 5.1, in order to validate the global estimate (2.21) 
we are left to show 

(5.13) \\ut(t a )\\ L 2 {Bl) < for all t > T Q . 

We will next show this is indeed the case. 
Lemma 5.3. Let u(x,t) be as in Theorem 1.6, then there exists Tq > such that 

(5.14) IM*o)IUa(Bi) < V^o for all t > T . 

Proof. Let To > 1 be as in (5.11) so that we have || lit (To) || £2 < ^/Eq. Since 
u(x,t) e W 1 ' 2 nC°°(B 1 x [l,oo),A0 and by the continuity of J B \u t (t)\ 2 in t, there 
exists 5 = S(Tq, £q) > such that for any to G [To, To + S] we have 

(5.15) lk(*o)|U=( Bl) < V^o~- 

Therefore by our previous arguments (especially Theorem 2.6, Lemma 5.1 and 
(5.12) with T) replaced by to), Lemma 2.5 applies to any subinterval of [To, To + 5] 
and yields 

(5.16) / \u t (t 2 )\ 2 < [ \ut(h)\ 2 for any t u t 2 G [T , T + 5}. 
JBx Jb 1 

This shows, instead of (5.15), for any to G [To, To + 8] we have 

(5.17) \\u t {to)\\ L 2 {Bl) < \\u t (T )\\ LHBl) < Vio". 

We can then continue and iterate this process beyond Tq + S and we see that 
Ib 1 l M *W| 2 i s indeed non-increasing along the flow after T). □ 

This completes the validation of (2.21) and therefore the assumptions in Lemmas 
2.4 and 2.5 in view of Theorem 2.6, finishing the proof of our main Theorem 1.6 as 
shown in Section 2. 
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